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Abstract 

Cardiac aging is characterized by diastolic dysfunction of the left ventricle (LV), which is due in part to increased LV 
wall stiffness. In the diastolic phase, myocytes are relaxed and extracellular matrix (ECM) is a critical determinant to 
the changes of LV wall stiffness. To evaluate the effects of ECM composition on cardiac aging, we developed a 
mathematical model to predict LV dimension and wall stiffness changes in aging mice by integrating mechanical 
laws and our experimental results. We measured LV dimension, wall thickness, LV mass, and collagen content for 
wild type (WT) C57/BL6J mice of ages ranging from 7.3 months to those of 34.0 months. The model was 
established using the thick wall theory and stretch-induced tissue growth to an isotropic and homogeneous elastic 
composite with mixed constituents. The initial conditions of the simulation were set based on the data from the 
young mice. Matlab simulations of this mathematical model demonstrated that the model captured the major 
features of LV remodeling with age and closely approximated experimental results. Specifically, the temporal 
progression of the LV interior and exterior dimensions demonstrated the same trend and order-of-magnitude 
change as our experimental results. In conclusion, we present here a validated mathematical model of cardiac 
aging that applies the thick-wall theory and stretch-induced tissue growth to LV remodeling with age. 



Introduction 

Over 70% of 50 million Americans over 60 years of age 
have cardiovascular disease (CVD) [1,2]. CVD preva- 
lence increases with age, and outcomes in older patients 
with acute coronary syndromes are poor [1,3]. It is 
important to note, however, that aging itself, even in the 
absence of CVD, alters LV structure and function and 
impairs the ability of the LV to respond to stress and 
injury. Thus, understanding the mechanisms of cardiac 
aging has significant clinical relevance. 

While most studies focus on the myocyte contribution 
to cardiac systolic function, indices for systolic function, 
such as ejection fraction, systolic velocities, and systolic 
isovolumic acceleration rate, have been shown to have 
little relation with age in both clinical and animal stu- 
dies [4,5]. In contrast, echocardiographic indices of 
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diastolic function including early (E) and late (Atrial- A) 
diastolic peak filling velocities and the E and A velocity 
ratio (E/A ratio), mitral deceleration time (the time 
from the peak to the end of the Doppler E-wave), isovo- 
lumic relaxation time (the time between the closing arti- 
fact of the aortic valve), and the earliest detection of 
trans-mitral blood flow have been demonstrated to 
decline with age in clinical studies [6-8]. With physiolo- 
gical aging, the LV undergoes monotonic structural 
changes that include increased wall thickness, chamber 
diameter, and mass [6-8]. 

We have shown diastolic dysfunction at the organ 
level in mice during cardiac aging [9], and diastolic dys- 
function is caused by increased myocardial stiffness at 
the tissue level [10]. The myocardium is composed of 
myocytes (muscle) surrounded by the ECM environ- 
ment. Accordingly, myocardial stiffness is determined by 
the volume ratio and the combined mechanical property 
of the myocytes and ECM. In the diastolic phase, myo- 
cytes are relaxed and ECM is a critical determinant to 
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the change of LV wall stiffness. About 90% of cardiac 
ECM composition in the young LV is collagen I and III. 
We have shown that collagen content in senescent 
mouse hearts doubles compared to young hearts [9,11]. 
Since collagen has a magnitude increase in stiffness over 
myocytes, the age-related increase of collagen content in 
the LV shifts myocardial mechanical properties from a 
myocyte-based stiffness to one influenced by collagen- 
based stiffness [12]. Therefore, the goal of this study 
was to evaluate the effects of ECM composition on LV 
remodeling with aging using a mathematical model 
developed by integrating cardiac mechanics and our 
experimental results in mice. 

Different LV wall stress models, such as the Laplace 
law based thin-wall models, thick-wall shell models, and 
finite element models, have been established to describe 
LV mechanics and compute stress [13,14]. Most current 
thin-wall and thick-wall models of LV remodeling were 
established assuming idealistic spherical, spheroidal, or 
ellipsoidal geometries. While finite element models 
allow some flexibility in the LV geometry, they require 
high computational power. Currently, these models 
focus on the stress calculation taken from a particular 
snapshot of the cardiac cycle or for the entire cycle. 
However, the interplay between LV stress and strain 
and between remodeling and geometric evolution on the 
life-time scale has not been established. 

Phenomenological models on geometric remodeling 
with aging have been developed to apply on arteries 
under hypertensive conditions [15,16]. The arterial wall 
is considered as a thick-wall vessel, and the remodeling 
equation on either inner or outer artery surface has 
been postulated as a function of strain and stress at the 
corresponding location. While these models do not con- 
sider the intrinsic relationships between phenomenologi- 
cal assumptions, they provide a possible methodology to 
model the temporal progression of tissue remodeling 
with aging. 

There are very few computational models available to 
study LV geometric adaptation with aging. Recently, we 
established a computational model of LV aging that 
incorporates Laplace law based stress model [11]. This 
model captures the overall trend of LV radius change 
with age. However, this thin-wall model assumes a con- 
stant thickness and no stress variation across the wall, 
LV thickening with age was not addressed. Here, we 
improve this model by adopting a more sophisticated 
thick-wall theory and stretch-induced tissue growth the- 
ory for the model used in this study. The novelty of this 
model lies in the integration of both computational and 
experimental approaches. The wall radii remodeling 
model established in this study was subject to the tem- 
poral function of the total mass which was measured in 



our experiments. Additionally, the predictions of LV 
geometry and end diastolic pressure-volume relation 
from the mathematical model were compared with our 
experimental measurements to validate predictions of 
the mathematical model. 

Methods 

With physiological aging, the LV undergoes tonic struc- 
tural and functional changes including increased mean 
LV wall thickness, chamber diameter, mass, concentric 
remodeling, and decline in LV diastolic function in 
humans [6-8]. We assumed mice age-related remodeling 
follow similar rules. Our model was established based 
on the following 2 assumptions. 

1) The LV was modeled as a concentric spherical shell 
with interior and exterior radii J?, and R 0 , respectively, at 
the no-load state. Due to the transmural pressure differ- 
ence P on the inner and outer surfaces, the deformed 
inner and outer radii were r; and r 0 , respectively. The 
ventricular wall was modeled as an isotropic and linear 
elastic material of two constituents, cardiac myocytes 
and ECM which is represented by collagen. 

(2) Internal and external radii, ie., R t and R ol as well as 
wall thickness R a - R h all have non-decreasing trend 
with age. 

Thick-wall LV model 

Based on assumption 1), deformation of the LV wall 
obeys Hook's law, which can be expressed in tensor 
form: e,j- = ^r(?ij — jjSijGkfo where L is strain tensor, a is 
stress tensor, v is Poisson's ratio, E is material's Young's 
modulus, 8jj is Kronecker delta, i, j, k e {R, 6, <f>} are the 
spatial indices for radial, meridional, and circumferential 
directions. The strain in the radial direction is: 

3r 1 

£rr = — ~ 1 = -[orr - v{0 00 + a^)]. (1) 

The stresses for elastic, isotropic, homogeneous, and 
spherical thick-wall contains the following analytical 
form[17]: 

PR 3 R 3 

PR 3 1 R 3 

aoe{R) = a^{R) = -^^{1 + -^), (3) 

in which, P denotes pressure difference applied to the 
inner and exterior walls, and R e [R h R 0 ]. Applying 
equations (2) and (3) to equation (1), we had 

r'^iWrT 111 - 2 "- 11 "*' 1 * 1 
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Integrate both sides of Eq. (4) from /?, to R ot we 
obtained 



r 0 -n = (l - 



E (R„/R,f 



K1 + v) Wf^ _ {1 _ 2vmK _ Ri) , (5) 



It's concluded that Poisson's ratios of all stable isotro- 
pic materials falls in the region (-1, 0.5) and only auxetic 
materials with honeycomb structures and networks have 
been found to have negative values [18]. Apparently, 
heart tissue is not auxetic. Therefore, for LV, its Pois- 
son's ratios falls in the region 0 <v < 0.5. As a conse- 
quence, it is apparent that the loaded wall thickness, i.e. 
r a - r t is smaller than the zero-load thickness, i.e. R a - R t . 
The LV transmural pressure P has been shown to alter 
with aging. We adopted the temporal profile of pressure 
difference reported in the literature [11,19]. 

The Young's modulus E increases with aging due to 
changes in the LV constituents. We previously reported 
that the fractions of the two constituents, i.e. myocytes 
and collagen changes with aging.[9,ll] In this study, we 
adopted a simplified version of the linear mixture theory 
of composite material, Young's modulus obeys the fol- 
lowing equation [20,21]: 



E{t)=E c v c {t)+E m [\-v c {t)], 



(6) 



in which E c and E m are two constants representing 
Young's moduli of collagen and muscle, respectively. 
Function v c (t) is the volume fraction of collagen, which 
can be expressed as: 



Vc{t) 



Vc{t) 



M c {t)/ Pc 



V c {t) + V m {t) M c {t)/Pc+M m {t)/ Pn 



(7) 



In equation (7), V c (t)/V m (t) is the volume occupied by 
collagen/muscle at time t, and M c (t)/M m (t) is the mass 
of collagen/muscle at time t. Parameters p c and p m are 
mass densities of collagen and myocytes. The density of 
collagen is an intrinsic property of the material and 
therefore does not change with time unless the collagen 
molecular structure undergoes significant changes. The 
density of collagen is adopted as p c = 1.70 g/ml [22]. 
For muscle, due to the fact that approximately 80% of 
the mass is water as has been measured in experiments 
for multiple mammalian muscle (either cardiac or non- 
cardiac), the density is close to that of water, with varia- 
tion contributed by tissue and tissue solids. We adopt 
myocyte density as p m = 1.06 g/ml [23,24]. 

Due to the fact that both collagen and muscle mass 
evolve with aging, the volume of the wall changed 
accordingly. The mass evolution enforces the following 
equations: 



M c (t) = Pcv c {t) \ 

jRi(t) 



4nR 2 dR, 



(8) 



M m {t) = p m [l - u c {t)] / 4nR 2 dR. 



(9) 



Take the derivative of equations (8) and (9) with 
respect to time, we obtained 

m c {i)Ipc = 4wv c (t)[R 2 „Wom - Rf (t)Ri(t)] + ^-[R 3 0 m - RfMRM. (10) 



M«(t)/ft» = - u c (t)][«?(t)R»W - Rf(t)«i(01 - x |R » (t) - R i Wl^W- (11) 

Add up Eqs. (10-11) and divide both sides by An, the 
following equation about free radii was obtained: 

MAt) + &-M m (i) , , 

R 2 0 Wo(t) - r 2 mm = — . (12) 



4irp c 



Wall incompressibility yields: 
r 3 0 (t)-rf(t)=R 3 0 {t)-Rl(t). 



(13) 



Since we assumed that the material is isotropic and 
homogeneous in our model (assumption 3), the spatial 
remodeling can only occur in the radial direction. As a 
result, the spherical geometry is always maintained. 
Therefore, only radii, e.g., the free radii i?, and R„ alter 
with aging. 

Stretch-induced tissue growth model 

The dynamic equation that governs the temporal beha- 
vior of the exterior radius in the free state was defined 
as <x(t) := which is the ratio of free exterior radius 

to its initial value. Further, since the deformed wall 
thickness is smaller than that in the free state, volume 
conservation between free state and deformed state of 
LV wall guaranties r a >R a and r, > R t . Additionally, the 
ratio between deformed and free exterior radii 
X[t) := jpjjj is always larger than 1 (X(t) > 1) since r a 

>R a . We proposed that the rate of change of free exter- 
ior radius had a power law relation kernel with the ratio 
X(t) and the curve followed an overall exponential decay, 

d(t) =f[t) = r R l [X D {t) - \\exp[-{t- t 0 )/z f ], (14) 

in which, / (t) is the function that characterized the 
rate of change for the outer radius, parameter x R is the 
remodeling characteristic time for R a (t), variable 

A D (t) - 1 = roD( ^uf t f ft) is the ratio of the change 

between deformed and free radii raised to the power D 
on the exterior surface, parameter Xf is the decay con- 
stant for function /which guarantees a bounded outer 
dimension of the LV. The parameters t R , D and tf are 
constant parameters to be determined. 

Equations (5-7) and (12-14), together with the defini- 
tions of a(t) and X(t), formed a mathematical model to 
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characterize the temporal progression of the LV geome- 
try by applying thick-wall theory and stretch-induced 
tissue growth. Given the temporal functions M c (t), M m 
(?) and P(t) with aging, we solved the mathematical 
model numerically. 

Input functions and parameters 

M c (t), M m {t), and P(t) are temporal functions established 
by fitting experimental data. We measured the LV mass 
of 148 wild type mice at different ages by necropsy. We 
adopted a smooth monotonic function that saturates with 
time to fit the experimental data for LV mass. The total 
LV mass function used in our simulation is M(t) = 76.88 
+ 8.92 {1 - exp[- (t - 8.18)/r iV } (solid line in Figure 1), 
the mass unit is milligram, the unit of time t is months, 
and r LV = 6.08 months. One of our recent studies showed 
that the collagen content (volume fraction) in wild type 
C57/BL6J mice is doubled in senescent mice (-1%) com- 
pared to young mice (-0.5%). Our previous studies also 
demonstrated the same trend of collagen content increase 
[11]. In our model, we chose 1% collagen content for 
young at 7.5 months, and approaching 2% for old mice 
after 30 months. By adopting the converted collagen con- 
tent, LV mass function M(t) above, and corresponding 
densities (see Table 1), the temporal curve for collagen 



Table 1 Constant parameters and initial values adopted 
in numerical calculation and the source of the data. 



Parameters (unit) Values References 



Initial value R,{t = 75mons) (mm) 


2.0 


Our experiments 
on over 148 mice 


Muscle density p m (gem" 3 ) 


1.06 


[23,24] 


Collagen density p c (gem -3 ) 


1.70 


[22] 


Muscle elastic modulus E m (Pa) 


5.0x1 0 4 


[28] 


Collagen elastic modulus E c (Pa) 


3.2x1 0 7 


[29,30] 


Poisson' ratio v 


0.47 


[31] 


Time step h (hrs) 


3 





mass was M c (t) = 1.21 + 1.50 { 1 - exp[- (t - 7.50) jx c \ 
with x c = 7.23 months and other units are the same to 
relevant ones in M (£). Since our model assumed that the 
LV consisted of two constituents (assumption 4), muscle 
mass were calculated by subtracting collagen mass from 
the total LV mass as M m (t) = M(t) — M c (t). Figures 2 
showed the temporal profiles of the mass of the total LV, 
collagen, and myocytes, as well as collagen volume frac- 
tion with age. Based on our computational prediction, 
total LV mass underwent a 12.8% increase from 7.5 to 30 
months. The mass growth for collagen was about 121% 
in parallel with a twofold increase in volume fraction as 
shown in Figure 2D due to a higher density in relative to 
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Figure 2 Computational results for mass functions for LV, collagen, and muscle, as well as volume fraction for collagen These 
functions were used as inputs to our mathematical model. (A) The total LV mass fitted with data generated in autopsy experiments; (B) Collagen 
temporal profile fitted with known trend from experiments; (C) Temporal profile of muscle mass calculated as the difference between total and 
collagen mass; (D) Volume fraction of collagen, the values for young (6-9 months) and senescent (>25 months) were obtained from published 
results of mouse. 



muscle. The myocytes mass grew about 11.1% in senes- 
cent mice compared to young group. 

Time-dependent pressure difference between the inter- 
ior and exterior surfaces of the LV was another input 
function. We collected the pressure measurement in the 
literature and derived its temporal progression as 
P(t) = 4.5 + 3.0 {l - exp[- (t- 2. 5)/r p ]} with unit 
mmHg (Figure 3A), where r p = 14.0 months [11,19]. The 
pressure ranged from 4.5 mmHg at 2.5 months to 7.1 
mmHg at 30 months. We also showed the calculated 
Young's modulus of LV wall by assuming linear and 
homogeneously mixed composite made of muscle and col- 
lagen (Figure 3B). The modulus increased by 82% in senes- 
cent mice compared to young mice. Though there is only 
an increase from 1% to 2% increase of collagen volume 
ratio, it gives rise to significant enlargement in the modu- 
lus of the composite because the elastic modulus of col- 
lagen is much larger compared that of the muscle. 

Further, from equations (2-3), the pressure difference P 
is the only external load that deforms the LV wall. With 
retreat of this pressure, or establishment of balance 
between inner and outer surfaces, any stress or strain 



disappears. Now if we compare the elastic modulus and 
the pressure difference, the highest pressure, i.e. 7.1 
mmHg, is equivalent to 944 Pa, less than 2 orders of mag- 
nitude lower the magnitude of Young's modulus. There- 
fore, without detailed computation, we can see that there 
should exist none noteworthy deformation (< 1%) and the 
strained dimensions should be extremely close to the free 
ones. This was proved by detailed computation shown 
later. 



Discretization of equations 

For convenience in the writing of equations, we defined 



c{t) 
c{t) 



pm 



m r°(o 3 

p(0 i 



I 2 



_ [{1 + v) whm 



(l-2u)])[R 0 (f)-Ri(t)], 
(l-2v)]}[R 0 (t)-R((t)l 



in which R°(t) = R 0 {t)/Ri{t), and d{t) = R 3 0 {t) - i?f (t)- 
Now Equations (5), (12), and (13) can be rewritten as: 



R 1 0 W 0 {t)-Rf{t)R i {t) = btt), 
r 0 {t) - r,{t) = c{t), 



(15) 



(16) 
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Time (Mon) Time (Mon) 



Figure 3 Temporal profile of pressure difference of LV wall and predicted Young's modulus of the myocardium (A) Time-dependent 
pressure difference across the LV. The curved is derived by fitting experimental data; (B) Elastic modulus of LV, which was calculated as isotropic 
and homogeneous composite made up of muscle and collagen. The modulus was calculated by taking account of temporal profiles of muscle 
and collagen in LV. 



r 3 0 {t)-rf{t)=d{t). 



(17) 



Given R t (t), R Q {t), P{t) and, E{t) variables c{t) and d(t) 
could be determined and r/f) and r B (t) can be solved as: 



r 0 {t) 



-c(t) + J[4d(t)/c{t) - c 2 (t)]/3 



c(t) + V[4d(0Mt)-c 2 (t)]/3 



(18) 



(19) 



The radii r, (t) and r 0 (t) were real positive numbers 
subject to the condition d{t) >c 3 (t). 

We studied the variation of free and deformed radii in 
the time domain [t 0 , t 0 + T\ and evenly divided this region 
T by a total of N small constant intervals of h. Equivalendy, 
the domain could be rewritten as [t 0 , t 0 + Nh]. The follow- 
ing is the procedure to calculate time- dependent radii. 

1. Read initial values Ri(t 0 ), R o (t 0 ), calculate r ; (£ 0 ) and 
r o (t 0 ) and A(t 0 ) = jf^p note that a(t 0 ) = 1. Iterate steps 
2-5 until the final step N. 



2. According to discretized equation (14) using Euler 
method, perform the following compute: 

X D (n) — 1 

a(n + 1) = a(n) + exp(—nh/tf)h, (20) 



R 0 {n+ 1) = a{n + l)R o {t 0 ). 



(21) 



3. Utilize discretized version of Eq. (15) taking 
account of Eq. (14) to calculate the R t at the next time 
step. 



Rcfn + 1) = Rifnl + l^44l 2 |A D (n) - 1 



(n) + l W) |2|iD(n) ~ ^(-"Vv^M - ^h. (22) 

4. Update values for c{n) and d{n) using their defini- 
tions by substituting the newly derived R t and R a . 



d{n + 1) =R 3 0 (n+ 1)-R?(n+ 1). 



(24) 
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5. Compute new values for r b r a and A, 

-c{n + 1) + y/\Ad\n + l)/c(n + 1) - c 2 (n + l)]/3 ^5) 



r,(n + 1) : 



r„(n+ 1) 



c{n + 1) + ^/[4d(n + l)/c(n + 1) - c 2 {n + l)]/3 



X(n+ 1) = 



r 0 (n+ 1) 
R 0 (n+1)' 



(26) 



(27) 



The numerical calculations were performed using 
in-house Matlab codes. In all simulations, parameter 
t R = 0.22 month, D = 2.5 and tj = 5.6 month, initial 
values used in our model were listed in Table-1. The 
initial radii were determined by minimizing the error 
from experimental measurements considering both 
end diastolic dimension/diameter (EDD) and wall 
thickness. 

Results 

Changes of EDD, wall thickness and stress 

The numerical simulation on LV dimensional changes 
including EDD and wall thickness were shown in Fig- 
ures 4. The interior radius increased from 2.00 mm to 
about 2.05 mm (see Figure 4B), which is consistent with 
the 2.60% increase in our mice experiments for EDD 
from young to senescence. None of the percentage dif- 
ference between strainless and deformed radii exceeds 
0.1% for all time of concern, indicating a proper adop- 
tion of small deformation Hook's law (see Figure 4C). 
The wall thickness started from 0.92 mm, ended with 
0.98 mm and had an increase of 0.053 mm, or a 5.72% 
ascent (Figure 4D), which is of the same magnitude as 
the 13.20% increase in our experiments. 

The temporal profiles of the mean o RR and mean Ogg 
(or a^,) across the wall were shown in Figures 4E and 
4F calculating by 



1 f R ° PR 3 R 3 



RijRo + lRi) 
~ 2{R2 + R 0 Ri + Rj)' 



(28) 



< (700 > = < 



R 0 

1 f PR 3 R 3 

R, (29) 
= Rj{R 2 0 + RpRj + 4Rf) 
4{R 3 -R 3 ) 

As is seen in Eq. (2), a RR < 0 and its contribution is a 
force toward inner direction. The magnitude of all these 
averaged stresses went up with time, which was a 



natural result of enhanced pressure P(t) (Figure 3A). We 
also showed that the magnitude of the circumferential 
or meridional stress was higher than that the radial 
(comparing Figures 4E and 4F). 

Pressure-volume relationship with aging 

End diastolic portion of pressure-volume loop in a car- 
diac cycle provides important information of passive fill- 
ing for LV and therefore describes passive properties of 
myocardium. For example, the slope of end diastolic 
pressure-volume relationship (EDPVR) gives the mea- 
sure of LV stiffness. We plotted end diastolic pressure- 
volume curve with age for C57/BL6J mice (see Figure 
5), in which EDV = (4/3)7rr 3 is the volume of the cavity. 
The arrow indicates the aging direction. The ratio of 
changes between pressure and end diastolic volume 
goes up with age, together with increased pressure and 
volume with age. It is worth noting that the end diasto- 
lic P-V relation shown in Figure 5 is a temporal progres- 
sion with age, which is different from the typical one 
cardiac cycle P-V loop. 

Analysis of constraints of remodeling equation 

We also analyzed the candidates of the tissue growth 
function f(t) based on assumption 2). From Eq. (12), we 
had: R 2 (t)R;(t) = R 2 (t)R 0 (t) - b{t). So R;(t) > 0 is 
equivalent as R 0 (t) > £>(t)R~ 2 (t). Apply dynamic equa- 
tion, i.e. equation (14), R„[t) = R o {t 0 )a(t), the inequality 



became: f(t) 



b(t) 



■pr. On the other hand, wall thick- 



Ro(t 0 )R 2 „(ty 

ness growth also gave R„(t) > Ri(t)- Since 



Rf(0 



Rf(0 



b(t) 



Ro(t) < 2 v ' 2 > substituting R a with a and using equa- 

K„ — R: 



tion for 

m< 



a, 

HO 



the 



R 0 {t 0 )[Rl{t) - Rf{t)] 



inequality 
, which led to 



Hi) 



Ro(t 0 )R o 2 (0 



<m< 



Hi) 



R o (t o )[R 2 (0-R 2 (t)]' 



(30) 



where b{i) was a universal input function in our sys- 
tem and did not depend on the selection of parameters. 
R o {t 0 ), i.e. initial value of R ol was a constant. However, 
R 2 (t) and R 2 (t) - R 2 (t) could only be determined after 
the simultaneous equations being solved. Therefore, 
these two terms were subject to the parameters 
adopted in the mathematical model. However, we 
could find universal upper and lower bounds for the 
growth rate function, J[t). As we know, R a (t) increases 

with aging. Hence, R j^j is a natural high upper bound 
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Figure 4 Temporal curves of LV dimensions and stresses derived by solving the simultaneous equations. The parameters in Eq (14) 
adopted are Tr = 0.22 mons, £) = 2.5 ar| d tf = 5.6 mons. (A) Temporal profile of the strainless radius on the external surface of LV; (B) The 
internal radius change with age; (C) Percentage change between deformed and free states. They are almost identical; (D) Wall thickness that has 
a 5.72% thickening from young to old mice; (E) Absolute value of negative spatial mean strain in radial direction; (F) Circumferential (or 
meridional) stress averaged across the wall. 



of fit). By defining S(t) = R%[t) - Rf{t), it is easy to see 
S(t) went up with time because both inner, outer radii 
and wall thickness enlarge. So S(t) was maximized at 
30 months which was the terminal simulation time 
span. To achieve the maximum surface area difference 
for two concentric spheres kept at a constant volume, 
the smaller the inner/outer radius, the bigger the value 
S is. Since both interior and exterior radii increase, the 
lowest possible radius at the internal surface was R t {t = 



30 mons) = Ri(t 0 ). Use the mass of each component at 
the final time and corresponding density, the final 
volume was determined and Parameter R a (t = 30 mon) 
was calculated, which is the exterior radius in month 
30 assuming the interior radius is the same as 
that in month 7.5, i.e. its initial value. Hence, 

■— — . ; : ^tt rr offers a lower 

R 0 (t 0 )[Rl(t = 30mons) -Rf(t = 30 mons)] 

upper bound. Therefore, if the following inequalities 
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up with age. 



are satisfied, the selection of the radius growth rate 
function fit) should give us proper remodeling trend, 



mo) 



6(0 



R D (t 0 )[R2(t = 30mons) - Rf(t = 30mons)\ 



(31) 



And the above inequalities do not rely on a concrete 
system. Both sides are known functions without para- 
meter dependence. Therefore, it is a universal condition 
for the selection of the type of rate function, fit) Though 
this is a stringent condition, functions fall in the pro- 
posed region will guarantee the stability of the system. 

As is shown in Figure 6, the green line with circle 
markers and red line with hexagram markers are the 
universal boundaries determined by equation (31). The 
blue curve with star markers and the black one with 
square markers are true upper and lower bounds respec- 
tively tor fit) calculated in equations (30) for the set of 
parameters listed above. 

In our simulations, parameters r R = 0.22 months, D = 
2.5 and if = 5.6 months were chosen by minimizing the 
deviations of computed EDD and wall thickness from 
measurements at senescent age, i.e. (/?, - Ri, exp ) 2 + [WT 
- WT exp ) 2 , in which WT means wall thickness. The 
growth function f(t) determined by these parameters, 
the dotted cyan curve, falls in the bounded region. The 
analysis of the boundary of the growth function 



provided a standard procedure for us to determine the 
possible sets for parameters with or without simulation 
results, which should be able to be applied in the mod- 
eling of age-dependent remodeling process of different 
species given experimental data, e.g. rat, dog or even 
human. In addition, if a growth function f(t) of this 
form is given, we can also use this boundary condition 
to predict whether LV geometry and wall thickness will 
monotonically increase or not. 

Discussion 

In this study, we established the first analytical mathe- 
matical model to quantify long-term LV geometric 
remodeling dynamics with aging by applying thick-wall 
theory and stretch-induced tissue growth postulate. This 
model addressed the temporal progression of cardiac 
aging and is an advancement from most current models 
that analyze the static equilibrium of LV remodeling 
and our previous thin-wall model of LV remodeling. In 
addition, we analytically determined the boundary for 
the tissue growth rate to guarantee the stability of the 
remodeling. This boundary might provide a reference 
for experimental measurement to examine the remodel- 
ing outcomes. Further, the parameters in the mathema- 
tical model were determined by using our experimental 
results from over 140 C57/BL6 mice. Some parameters 
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and initial conditions were selected by minimizing the 
error between our computational predictions and 
experimental measurements. In addition, the model was 
validated by comparing the predictions of internal and 
external radii and wall thickness in diastole to experi- 
mental results. The predicted LV geometry trends were 
consistent with the experimental results. Thus, this 
study is a real integration of computational and experi- 
mental approaches for model establishment and valida- 
tion using mice data. This approach and the proposed 
model can be applied to establish a cardiac aging model 
for human in the future. 

The P-V curve in Figure 5 shows that both pressure 
and end diastolic volume escalate with age. This has 
been observed in a LV P-V relationship experiment on 
rat aging study by Pacher and colleagues [25]. In their 
study, LV end diastolic volume, LV end diastolic pres- 
sure, and the slope of EDPVR were greater in senescent 
than in young rats. In another experimental aging study 
on C57BL/6 female mice [26], the pressure was augmen- 
ted in 16-month-old mice compared to 6-month-old 



group. In our experimental study, we observed increased 
LV volume with age. These experimental results con- 
firmed our computational predictions. 

Though our computational results follow the trend of 
LV geometry changes, the proposed model has some 
limitations. First, we assumed the LV to be a spherical 
shell. With this assumption, all directions grow at the 
same pace to maintain the spherical symmetry. Sphere 
is a special geometry that requires minimum linear 
dimension for a constant volume. However, LV does not 
have a spherical geometry in reality. In addition, we 
assumed that LV was composed of elastic, isotropic, and 
homogeneous materials and applied the linear mixture 
theory to calculate the Young's modulus of the myocar- 
dium. Though the assumptions simplify the mathemati- 
cal model and the computational analysis, such 
simplification might lead to prediction errors. From our 
computational results, the radii growth calculated are 
3.43% for external, 2.37% for internal radius, and 5.72% 
for wall thickness collected from C57/BL6 mice. Some 
of these values were small compared to our echo 
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experiment on C57/BL6 mice, in which EDD grows 
2.60% and wall thickness is up 13.20% from young to 
senescent. In addition, our previous experimental work 
on CB6F1 mice reported EDD and wall thickness have 
12.73 and 13.00% growth, respectively, in the LV 
between the groups of young and senescent [27]. This 
discrepancy may be originated from the simplification of 
the geometry we chose or the prediction of the myocar- 
dium stiffness. Future experiments will build on these 
limitations to improve our results by incorporating 
more realistic geometry and taking more advanced tech- 
niques to determine LV stiffness. 

Conclusions 

We have established the first mathematical model to 
study age-related temporal-spatial LV remodeling by 
adopting thick-wall theory and stretch-induced tissue 
growth theory. Ilnputs of the mathematical model were 
real experimental data including temporal profiles of LV 
mass, collagen content change, and pressure across LV, 
which were obtained by over 140 mice. The established 
model captured the major property of LV remodeling 
with age and yielded predicted results of LV geometry 
progressions in mice comparable to experimental results. 
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